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Abstract. We review our investigation of the spin dynamics for two electrons 

confined to a double quantum dot under the infiuence of the hyperfine interaction 

between the electron spins and the surrounding nuclei. Further we propose a scheme 

pH ' to narrow the distribution of difference in polarization between the two dots in order 

I ' to suppress hyperfine induced decoherence. 
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1 Introduction 



The fields of semiconductor physics and electronics have been successfully 

combined for many years. The invention of the transistor meant a revolu- 

■ O ■ tion for electronics and has led to significant development of semiconductor 

^ I physics and its industry. More recently, the use of the spin degree of freedom 

fj ■ of electrons, as well as the charge, has attracted great interest 1 . In addition 

to applications for spin electronics (spintronics) in conventional devices, for 
instance based on the giant magneto-resistance effect 2 and spin-polarized 
►^ I field-effect transistors 3 , there are applications that exploit the quantum 

fv^ ■ coherence of the spin. Since the electron spin is a two-level system, it is a 

lO i natural candidate for the realization of a quantum bit (qubit) T. A qubit 

Cn ■ is the basic unit of information in quantum computation. The confinement 

^i I of electrons in semiconductor structures like quantum dots allows for bet- 

vQ ■ ter control and isolation of the electron spin from its environment. Control 

f~^ I and isolation are important issues to consider for the design of a quantum 

computer. 

C^ I Formally, a quantum computation is performed through a set of transfor- 

mations, called gates 5 . A gate applies a unitary transformation U to a, set 
of qubits in a quantum state \'F). At the end of the calculation, a measure- 
ment is performed on the qubits (which are in the state \'P') = U |!Z')). There 
Q I are many ways to choose sets of universal quantum gates. These are sets of 

O ' gates fronr which any computation can be constructed, or at least approxi- 

mated as precisely as desired. Such a set allows one to perform any arbitrary 
. , ■ calculation without inventing a new gate each time. The inrplementation of 

rS ' a set of universal gates is therefore of crucial importance. It can be shown 

H : 

Cu * Presented as plenary talk at the annual DPG meeting 2006, Dresden (to appear 

in Advances in Solid State Physics vol. 46, 2006). 
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that it is possible to construct such a set with gates that act only on one or 
two qubits at a time |S]. 

The successful implementation of a quantum computer demands that 
some basic requirements be fulfilled. These are known as the DiVincenzo 
criteria |7] and can be summarized in the following: 

1. Information storage-the qubit: We need to find some quantum property 
of a scalable physical system in which to encode our bit of information, 
that lives long enough to enable us to perform computations. 

2. Initial state preparation: It should be possible to set the state of the qubits 
to before each new computation. 

3. Isolation: The quantum nature of the qubits should be tenable; this will 
require enough isolation of the qubit from the environment to reduce the 
effects of decoherence. 

4. Gate implementation: We need to be able to manipulate the states of in- 
dividual qubits with reasonable precision, as well as to induce interactions 
between them in a controlled way, so that the implementation of gates is 
possible. Also, the gate operation time Tg has to be much shorter than 
the decoherence time T2, so that T5/T2 ^ r, where r is the maximum 
tolerable error rate for quantum error correction schemes to be effective. 

5. Readout: It must be possible to measure the final state of our qubits once 
the computation is finished, to obtain the output of the computation. 

To construct quantum computers of practical use, we emphasize that the 
scalability of the device should not be overlooked. This means it should be 
possible to enlarge the device to contain many qubits, while still adhering 
to all requirements described above. It should be mentioned here that this 
represents a challenging issue in most of the physical setups proposed so far. 



2 Quantum-Dot Spin Qubit 

The requirement for scalability motivated the Loss-DiVincenzo proposal ^ 
for a solid-state quantum computer based on electron spin qubits. 

The qubits of the Loss-DiVincenzo quantum computer are formed from 
the two spin states (|T) j |i)) of a confined electron. The considerations dis- 
cussed in this proposal are generally applicable to electrons confined to any 
structure, such as atoms, molecules, etc., although the original proposal fo- 
cuses on electrons localized in quantum dots. These dots are typically gener- 
ated from a two-dimensional electron gas (2DEG), in which the electrons are 
strongly confined in the vertical direction. Lateral confinement is provided 
by electrostatic top gates, which push the electrons into small localized re- 
gions of the 2DEG (see Fig.^. Initialization of the quantum computer can 
be achieved by allowing all spins to reach their thermodynamic ground state 
at low temperature T in an applied magnetic field B (i.e., virtually all spins 
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Fig. 1. Two neighbouring electron spins confined to quantum dots, as in the Loss- 
DiVincenzo proposal. The lateral confinement is controlled by top gates. A time- 
dependent Heisenberg exchange coupling J(t) can be pulsed high by pushing the 
electron spins closer, generating an appreciable overlap between the neighbouring 
orbital wave functions. 



will be aligned if the condition jg^B^I ^ k^T is satisfied, with (/-factor g, 
Bohr magneton /ie, and Boltzmann's constant fee). Single-qubit operations 
can be performed, in principle, by changing the local effective Zeeman inter- 
action at each dot individually. To do this may require large magnetic field 
gradients |H1, 5- factor engineering [Hj, magnetic layers, the inclusion of nearby 
ferromagnetic dots ^, polarized nuclear spins, or optical schemes. 

In the Loss-DiVincenzo proposal, two-qubit operations are performed by 
pulsing the electrostatic barrier between neighboring spins. When the barrier 
is high, the spins are decoupled. When the inter-dot barrier is pulsed low, 
an appreciable overlap develops between the two electron wave functions, 
resulting in a non-zero Heisenberg exchange coupling J. The Hamiltonian 
describing this time-dependent process is given by 

H{t)^J{t)SL-SR. (1) 

This Hamiltonian induces a unitary evolution given by the operator U = 
Texp{—i J H{t)dt/h}, where T is the time-ordering operator. If the ex- 
change is pulsed on for a time Tg such that / J{t)dt/h = JqTs/Ti = n, the 
states of the two spins, with associated operators S_l and S_r, as shown in fig- 
ure Q] will be exchanged. This is the SWAP operation. Pulsing the exchange 
for the shorter time Ts/2 generates the "square-root of SWAP" operation, 
which can be used in conjunction with single-qubit operations to generate 
the controUed-NOT (quantum xor) gate 0] . The "square-root of SWAP" gate 
has recently been implemented in an experiment by Petta et al. |10| with a 



4 D. Klauser et al. 

switching time r^ = 180 ps. For scalability, and application of quantum error 
correction procedures in any quantum computing proposal, it is important to 
turn off inter-qubit interactions in the idle state. In the Loss-DiVincenzo pro- 
posal, this is achieved with exponential accuracy since the overlap of neigh- 
boring electron wave functions is exponentially suppressed with increasing 
separation. A detailed investigation of decoherence during gating due to a 
bosonic environment was performed in the original work of Loss and DiVin- 
cenzo. Since then, there have been many studies of leakage and decoherence 
in the context of the quantum-dot quantum computing proposal. 

In addition to the interaction-based gate operations introduced above, 
it has been shown recently |11[ [T^ that it is also possible to generate the 
controUed-NOT based on partial Bell state (parity) measurements. 

For both interaction-based and measurement-based quantum computa- 
tion with the quantum-dot spin qubit, decoherence due to the coupling of 
the qubit to its environment is a major obstacle. There are two important 
sources of decoherence in GaAs quantum dots: spin-orbit coupling (interac- 
tion between spin and charge fluctuations) and hyperfine coupling (interac- 
tion between the electron spin and nuclear spins). In the case of spin-orbit 
interaction alone it has been shown that the decoherence time T2 (which is 
the relevant timescale for quantum computing tasks) exceeds the relaxation 
time Ti and is given by T2 = 2Ti to leading order in the spin-orbit coupling 
|1H| . Since the Ti obtained in measurements ^1 El is on the order of ms, 
but the ensemble-averaged dephasing time T2* measured is ~ 10 ns, spin-orbit 
interaction is not limiting for the dephasing time T2*. The limiting source of 
decoherence is the hyperfine interaction ,10, . 



3 Hyperfine Interaction in Single and Double Dots 

The hyperfine interaction between the electron spin and the nuclear spins 
present in all III-V semiconductors ^Hl leads to the strongest decoherence 
effect [iniinillilllEniEIllllllSllMI- Experiments HOI 123 Ei EZI have 
yielded values for the free-induction spin dephasing time Tj* that are consis- 
tent with T* - ^/N/A - 10ns [1312111221 for N = 10^ ?i = 1, and A = gO/xeV 
in GaAs, where N is the number of nuclei within one quantum dot Bohr ra- 
dius and A characterizes the weighted average hyperfine coupling strength in 
GaAs [211 ■ This is to be contrasted with potential spin-echo envelope decay 
times, which may be much longer |28l 12911^(11 18 1| . With a two-qubit switching 
time of Ts ~ 180 ps ^01 this only allows ~ 10^ gate operations within T2, 
which falls short (by a factor of 10 to 10^) of current requirements for efficient 
quantum error correction |32) . 

There are several ways to overcome the problem of hyperfine-induced 
decoherence, of which measurement and thus projection of the nuclear spin 
state may be the most promising (23- Other methods include polarization 
|17l 1221 1281 ESI of the nuclear spins and spin echo techniques |10l 1281 IHTI) . 
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However, in order to extend the decay time by an order of magnitude through 
polarization of the nuclear spins, a polarization of above 99% is required 
123), but the best result so far reached is only ^--^60% in quantum dots |25) . 
With spin-echo techniques, gate operations still must be performed within the 
single-spin free-induction decay time, which requires faster gate operations. 
A projective measurement of the nuclear spin state leads to an extension of 
the free-induction decay time for the spin. This extension is only limited by 
the ability to do a strong measurement since the longitudinal nuclear spin in 
a quantum dot is expected to survive up to the spin diffusion time, which is 
on the order of seconds for nuclear spins surrounding donors in GaAs |21] . 

A detailed analysis of the spin dynamics for one electron in a single quan- 
tum dot can be found in Ref.^Hl- Here we concentrate on the case of two 
electrons in a double quantum dot. The spin S/ of an electron in quantum 
dot I = 1,2, interacts with the surrounding nuclear spins Ik via the Fermi 
contact hyperfine interaction: 

HM^Si-hi; h,=^4lfc; Ai ^ voA\^\rk)\\ (2) 

A; 

where vq is the volume of a crystal unit cell containing one nuclear spin. 
The effective Hamiltonian in the subspace of one electron on each dot is best 
written in terms of the sum and difference of electron spin and collective 
nuclear spin operators: S = Si -I- S2, i5S = Si — S2 and h = i(hi -|- h2), Sh = 
|(hi-h2): 

H.ffit) = ^zS' + h • S + 5h-(5S + ^S • S - J(i), (3) 

where e^ = g^sB is the Zeeman splitting induced by an applied magnetic 
field B = {0,0, B),B > 0. We assume that the Zeeman splitting is much 
larger than ((5h)i.ms and (hi)i.,„s, where (O)rms = {'4'i\ C^ li'i) is the root- 
mean-square expectation value of the operator O with respect to the nuclear 
spin state \ipi). Under these conditions the relevant spin Hamiltonian becomes 
block diagonal with blocks labeled by the total electron spin projection along 
the magnetic field S^. In the subspace of S*^ = the Hamiltonian can be 
written as (h ^ 1) ^1155] 

Ho{t)^^T''-^[2T'; J{t)^Jo+jcos{ujt), fl ^2{5h' + 5b'). (4) 

Here, (56^ is the inhomogeneity of an externally applied classical static mag- 
netic field with (56^ <C B. The Pauli matrices t", a = x,y,z are given in the 
basis of l-f) = |r- = 1) = lit) and |-> = \t' = -I) = |ti). 

The dynamics of the two-electron spin states depends strongly on the 
initial state of the nuclear spin system. We denote by \n) the eigenstates of 
6h' with Sh'^ \n) = Shf^ \n). If the initial state of the nuclear spin system 
is p/(0) = \n) {n\ and if we neglect spin-flip processes (as can be done for 
a large enough magnetic field B), then the initial spin state of the electron 
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does not decay. Thus, if it is possible to prepare the nuclear spin system in an 
eigenstate \n), hyperfine- induced decoherence could be overcome. In general, 
however, the initial state of the nuclear spin system is not an eigenstate |n) 
but a general mixture: 

pm^J2P^\i^\){i^]\; |^l)-5]a;|«), (5) 

i n 

where the a^ satisfy the normalization condition J2n l^nl^ = ^ ^^^ J2iPi — 1- 
We denote by pi{n) = J2iPiWn\'^ the diagonal elements of the nuclear spin 
density operator. 

For a large number of nuclear spins N ^ 1 which are in a superposition 
of (5/i^-eigenstates |n), pi{n) describes a continuous Gaussian distribution of 

Shf^ values, with mean 5h^ and variance cr^ = ((5/i^ — ShPj . In the limit of 
large N, the approach to a Gaussian distribution for a sufficiently randomized 
nuclear system is guaranteed by the central limit theorem |52]. We perform 
the continuum limit according to 

Yl /'^H/H -^ / dxpi{x)f{x); pi{x) = -j=- exp f ^^^° j ,(6) 

where x = Sh^ + Sb^, xq — 5h^ + 6b^ and Uq = x'^ — x'q. 

For the case of a static exchange interaction J{t) = Jq, the decay of 
the two-electron spin states in the S^ = subspace due to the Gaussian 
distribution of nuclear spin states may be calculated in several interesting 
limits I^ESl- Assuming the initial state of the two-electron system is Pe{0) = 
\+) {+\, the probability P+ to measure the |-|-) state as a function of time is 
given by 



2 s{xy 2s{xy 



with s{x) — vJ2 + 4a?^. In the limit co ^ 0, which corresponds to one fixed 
eigenvalue, there is no decay. However, for ctq > there is decay. For the 
regime |a:o| ^ ctq we have a Gaussian decay at short times with a decay time 
to ~ I/ctq: 

^i.(')-5 + f + (5-f)»p(-|)co,M, (8) 

3/2 

uJo^J.P+^xl to= /^" ; |a;o|>CTo, ^ < :7^- (9) 

V " 4|a;o|cro 2JVg 

Thus, decreasing ctq increases the coherence time Iq. Hence, the strategy to 
suppress hyperfine-induced decoherence is to narrow the Gaussian distribu- 
tion of nuclear spin eigenvalues through a measurement of the eigenvalue of 
6h^, i.e., of the difference in polarization between the two dots |23[I24| . It has 
also been proposed to measure the nuclear spin polarization using a phase 
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estimation method |36j . In the ideal case, phase estimation yields one bit of 
information about the nuclear spin system for each perfectly measured elec- 
tron. Optical methods have also been proposed [37] ■ The all-electrical method 
we propose here can be applied with current technology used in Refs. ;10..27) . 



4 Nuclear Spin State Narrowing 

The general idea behind state narrowing is that the evolution of the two- 
electron system is dependent on the collective nuclear spin state and thus 
knowing the evolution of the two-electron system determines the nuclear spin 
state. 

The eigenstates of the Hamiltonian Hq are product states: if the nuclear 
spin system is in an eigenstate \n) of Sh^ with 6h^ \n) = Sh'^ \n), we have 
H \tp) = Hn \ipe) (8) \n), where in i?„ the operator 5h^ has been replaced by 
Shf^ and \ipe) is the electron spin part of the wave function. Thus, in the 
Hamiltonian for the evolution of the initial two-electron system, the param- 
eter Shf^ is determined by the state of the nuclear spin system. Initializing 
the two-electron system to the |-|-) state, i.e., pe(0) — \+) (+| and performing 
a measurement in the |±) basis at time t„i yields for the distribution of nu- 
clear spin eigenvalues (which is the diagonal part of the nuclear spin density 
operator in the continuum limit) after the measurement |35j 



Pi 



^^'^+'-\x)^pj{x){l-LU^))^, (10) 



P. 



p^l'-'-\x) ^ pi{x)LUx)^, (11) 

where pi{x) is the initial Gaussian distribution of nuclear spin eigenvalues 
(see Eq. ©) and the probabilities P for measuring |±) are given by 

/•OC 

P+ = / dxpi{x){\~L^{x)), (12) 



P- - / dxpi{x)L^{x), (13) 

with 

1 (j74)^ 

"^-^"^ ^ 2(.-|)^ + (,74)- (^^) 

To obtain this result we have assumed that the measurement is performed 
with low time resolution [^ At ^ 1/j and that the parameters satisfy 
the requirements given in Eq. H15|l below. The distribution of nuclear spin 
eigenvalues after the measurement depends on the result of the measurement 
(whether |-f) or |— ) was measured) and on the driving frequency to of the 
oscillating exchange interaction J{t). In the case where the measurement out- 
come is |-|-), the initial distribution pi{x) is multiplied by 1 — Li^{x) which 
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causes a dip in pi{x) aX x — lo/2. However, in the case where the resuh 
of the measurement was |— ), the initial distribution pi{x) is muhiphed by 
L{x). The full-width at half-maximum (FWHM) of i^(x) is j/2, i.e., half the 
amplitude of the applied oscillation exchange interaction J{t). Thus, choos- 
ing j < ctq, Pj ' '" (a;) is dominated by the Lorentzian and therefore the 
FWHM of the initial nuclear spin distribution has been narrowed by a factor 
« j/4cto. The probability Pj to measure |— ) in the regime j <C cfq is given by 
Pj ~ (j/ScTo) exp((a;o — u!/2)'^/2(Tq) and the nuclear spin distribution after 
measuring |— ) is centered around lu/2. Thus, through such a measurement 
it is possible to choose the center of the nuclear spin distribution by choos- 
ing the driving frequency. However, the larger the difference xq — lo/2, the 
smaller is the probability to have measurement outcome |— ), which leads to 
narrowing. 

4.1 Experimental Recipe 

An experimental implementation of this scheme of course requires the ability 
to initialize to the state |-|-) and to read-out the states |±). This has recently 
been achieved in an experiment by Petta et al. QUI using adiabatic passage 
from the S^ — singlet. What needs to be achieved in addition is to apply an 
external magnetic field gradient Sb^ between the two dots in order to satisfy 
the requirements on the parameters of the system: 

-'^o < xq, j < a;o, ctq < xq, j < ctq. (15) 

Typical values for the parameters satisfying these requirements are: 1/cro = 
10 ns, l/j = 100 ns, lu ^ 2xo = 10^ - lO^'' Hz. 

The pulse-sequence for one measurement is shown in Fig. El The pa- 
rameter e describes the detuning between the singlet state with two elec- 
trons on the right dot and the singlet state with one electron on each 
dot: e = i?s(o,2) — -E-sCi,!)- First the system is set to the 5'(1,1) from the 
^(O, 2) state by going from large positive to negative detuning e (such that 
still J » \xo\) as described in Ref. JOI (rapid adiabatic passage through 
5(1, 1) — T+ resonance) . In the hmit of J ^ |a;o| and xo > 0, the ground 
state is |-|-) (for xq < 0, the ground state is |— ) and |±) thus need to be 
interchanged in the following description) and initialization to |+) is thus 
possible by adiabatic passage from 5(1, 1), i.e., by switching adiabatically 
to large negative detuning (such that J -^ xq). Then the oscillating signal is 
applied to J{t) for a time tm- Finally we adiabatically switch back to J ^ xq. 
With this the |-|-) state goes to the singlet S{1, 1), and the |— ) state goes to 
the S^ = triplet To(l, 1). Read-out of the singlet and triplet may then be 
achieved via switching to large positive detuning: the ^(l, 1) state goes over 
to the S{0, 2) while the ro(l, 1) does not since tunneling preserves spin. The 
number of electrons on the right dot can then be detected via a charge sensor 
(QPC). If the outcome of the measurement is |— ), we have achieved narrow- 



Quantum-Dot Spin Qubit and Hyperfine Interaction 




ii +) Kg(My 




Measure 





^x 


\pi(.x) 


/ 




\ 

\ 


/ 


pr 


,-2x„)(\\^ 


, , K 





(x - xa)/aa 



Fig. 2. In this figure the pulse-sequence for one measurement in the basis |±) 
is explained, (a) The level diagram for the two-electron spin states (sweeping 
e = -Es{i,i) — Es(o,2) with -Es(i,i) 4- Es(o,2) held constant). Inset: The splitting 
between the S(l, 1) and To state is given by the exchange interaction J between 
the two dots. For J ^ the j±) states become eigenstates. (b) The change of the 
detuning e during the course of the measurement is sketched: the position of the 
boxes corresponds to the value of e. After applying the oscillating signal the system 
is in either one of |±), which results in a different state when switching back to 
positive detuning. If the system ends up in the To state (which corresponds to mea- 
surement result |— )) narrowing has been achieved, otherwise the nuclear system 
must be rerandomized and the measurement repeated. 



ing. In the case where we have measured |+), the nuclear spin distribution is 
rerandomized and the measurement is repeated. 



4.2 Adaptive scheme 

If measurements at several different driving frequencies can be performed, 
a systematic narrowing of the distribution can be achieved by an adaptive 
scheme. Such an adaptive scheme is more intricate than the one described 
above, but allows one to narrow by more than a factor j /Agq. 

The results of Eqs. (|IOIII|l may be generalized to the case of M subsequent 
measurements at different driving frequencies uji under the assumption, that 
the nuclear spin system is static between subsequent measurements: 
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Fig. 3. In this figure we sliow a typical |39| sequence of the rescaled probability 
density ot eigenvalues tv(xj — pj (x)/max I pj (x) I tor the 



adaptive scheme. Here, pj 



{M,{a-},{i^,}) 



{x) is given in Eg. HIGH . We have x = Sh^ + SV 



j/(jQ — 1/10, a ~ '^ -I'^i J fl'iid in a)-c) the initial Gaussian distribution (with 
FWHM 2(To\/21n2 ~ 2ao) is plotted for reference, a) Up to M — 37 measurements 
the outcome is never |— ) and thus each measurement "burns a hole" into the dis- 
tribution where it previously had its maximum, b) In the 38th measurement the 
outcome is |— ) which leads to a narrowed distribution of nuclear spin eigenvalues 
(peak centered at « —0.25) with a FWHM that is reduced by a factor ~ j/Aao- 
c) Adapting the driving frequency uj to this peak, i.e., setting uj/2 = Xmax in sub- 
sequent measurements, leads to further narrowing every time |— ) is measured. In 
this example the final FWHM is « cro/lOO, i.e., the distribution has been narrowed 
by a factor ~ j/lOao- d) The probability P~ to measure |— ) jumps up after the 
38th measurement and after |— ) is measured several more times, this probability 
saturates close to 1/2. 



^(M,|. 
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(16) 



where Q({a~}, {uJi}) is the normahzation factor, a~ = 1 for measurement 
outcome |— ) and a~ — for measurement outcome |-|-) in the i* mea- 
surement with driving frequency cui. Further, {uji} = {lui, . . . ,ujm} and 
{a~} = {aj", . . . ,aj/}. As we have seen in the case of just one measure- 
ment, it is the measurement outcome |— ) that leads to narrowing. Thus, 
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before each measurement w^ is chosen to maximize the probabihty P^ to 
measure |— ). The reason that LOi must be adapted and that one should not 
keep measuring at the same driving frequency is that the measurement out- 
come |-|-) causes a dip in pi{x) at the position where L^.{x) has its peak and 
since P~. is the overlap of p/(a;) and L^^{x), this causes Pj". to diminish with 
each measurement. 

To see what is a typical measurement history for such an adaptive scheme 
we have performed simulations. The results for a typical [231 sequence of 
measurements is shown in Fig. O (for another sequence see Fig. 2 in Ref. [35! ) . 



5 Conclusions 

We have reviewed our scheme that uses pseudospin measurements in the 
S^ = Q subspace of two electron spin states in a double quantum dot to 
narrow the distribution of difference in nuclear polarization between the two 
dots. A successful experimental implementation of this scheme would allow 
to suppress hyperfine-induced decoherence and thus to reach the coherence 
times required for efficient quantum error correction. 
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